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We investigate the influence of order parameter fluctuations on the transition between normal and 
mixed superconducting states at low temperatures. We show that in case of clean quasi-two- 
dimensional superconductors the transition can be described by the functional of the Ginzburg- 
Landau type. We consider the large— N generalization of this functional and using the lowest 
Landau level approximation we get the large— N equations which describe the phase transition. In 
case of physical dimensionality we found that the transition is of the first order. The fluctuations 
significantly affect the temperature dependence of the upper critical field. 



I. INTRODUCTION . 

It is well known that the magnetic field penetrates 
type—// superconductors through flux lines which form 
the Abrikosov latticefil. The theory of this mixed state 
was first developed by Abrikosov for temperatures close 
to T c El, and then it was extended to all temperatureso. 
These theories ignore the fluctuations of the order param- 
eter. It is a very good approximation for the conventional 
superconductors because in this case the order parame- 
ter fluctuations are important only in a very small region 
near the phase transition line. This happens because the 
Ginsburg numbers of the conventional superconductors 
are very small. 

However for high— T c superconductors there are exper- 
iments which cannot be explained by the usual mean field 
theory. The upper critical field H C 2 at low temperatures 
significantly increases as temperature decreases instead 
of being approximately constant as follows from the mean 
field theorycl. The Ginzburg numbers of high— T c super- 
conductors are not very small, therefore one can suggest 
that this unusual behavior is due to the order parame- 
ter fluctuations. Also, it is important to understand the 
type of the phase transition (first order or second order) , 
because in mean field approximation this transition is al- 
ways of the second order, but the fluctuations can induce 
the first order phase transition. This happens, for exam- 
ple, in the model described in Ref U arid it was suggested 
to happen in the Abrikosov transitionB. 

In the following we will consider the influence of the 
order parameter fluctuations on the transition between 
normal and mixed superconducting states in clean su- 
perconductors. We argue that in case of pure quasi-two- 
dimensional superconductors, when the magnetic field is 
applied along the low-conducting direction, even at low 
temperatures, this transition can be investigated by the 
effective functional of the Ginzburg-Landau (GL) type 
which contains imaginary time, because quantum fluc- 
tuations become important at low temperatures. Un- 
fortunately, even having the effective functional one can 
hardly calculate the free energy exactly, which is typical 
for the critical phenomena problems. Therefore we mod- 



ify the functional introducing rt— index to the fluctuating 
field and considering the large— TV limit. One should be 
careful in introducing of n— index into this functional, 
because doing that in not a proper way one can get a 
model which does, not have a solution in the form of the 
Abrikosov latticed. Under the approximation when the 
mean field solution and fluctuations belong to the lowest 
Landau level (which is valid near the phase transition) we 
get the equations which describe the phase transition. 

We show that in high enough dimensions (either T = 
0, d > 4 or T 7^ 0, d > 6) if the coupling constant is 
not too large the transition is of the second order. In 
this case the corrections from fluctuations do not mod- 
ify mean field results essentially. In the physical case 
(d=3) the phase transition is always of the first order and 
the fluctuations significantly affect the phase transition 
line: the upper critical field increases as temperature de- 
creases, and the curvature of this dependence is negative 
when the temperature is not too low (see later, Figs. 4, 5). 

The paper is organized as follows: In Sec. || we de- 
scribe the model. In Sec. Ill we derive the large- N equa- 
tions describing the phase transition. In Sec. |y| we 
analyze the phase transition from the side of the normal 
state region. In Sec. [v] we simplify the large- TV equations 
making the lowest Landau level approximation. In Sec. 
VI we solve the l arge — A equations in case of high di- 
mensions. In Sec. VII we solve the large- A equations in 
the physical case. In Sec. VIII we discuss the spectrum 



of the order parameter fluctuations. And we summarize 
our results in Sec. [X. 



II. DESCRIPTION OF THE MODEL. 

To analyze the transition to the superconducting state 
one can usually use GL functional. This functional is 
based on the expansion in the order parameter, so it is 
valid when the ratio of the order parameter and temper- 
ature A/T is small. In case when an external magnetic 
field close to H c2 is applied, one can still use GL func- 
tional if A/T <C 1, but at low temperatures the expan- 
sion in this parameter becomes not possible. From the 
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other hand the magnetic field produces a depairing ef- 
fect, as temperature does, so one can try to expand in 
A because it is small compared with the magnetic field 
taken in the proper units. This was done in RefE and it 
was shown that the coefficient in |</>| 4 term has a logarith- 
mic singularity In i , which means that the expansion in 
A is not possible at zero temperature. The appearance 
of this singularity can be easily understood physically: 
Indeed, if we consider a pair of electrons which move 
exactly parallel to the magnetic field, then the magnetic 
field does not affect them in the quasiclassical approxima- 
tion. Therefore this singularity comes from the electrons 
moving parallel to the magnetic field or close to this di- 
rection, because the magnetic field does not produce a 
depairing effect for these electrons. The direct calcula- 
tion of the coefficient of |0| 4 term shows, that indeed, 
the singularity comes from the momenta parallel to the 
magnetic field. 

There is a way to avoid this problem if one considers 
a quasi- two-dimensional superconductor. Indeed, if the 
magnetic field is applied along the low-conducting direc- 
tion, then due to the quasi-two-dimensional band struc- 
ture there are no momenta parallel to the magnetic field. 
In this case it is possible to expand in A and get a func- 
tional of GL type. Unfortunately, the functional which 
arises does not match exactly the form of the usual GL 
functional. The difference is that the |0| 4 term becomes 
a nonlocal functional of (j> fields with a range of "inter- 
action" of the order of the magnetic length!! We hope 
that this difference is not crucial, and therefore we con- 
sider the model with a local |</>| 4 term. The Lagrangian 
density of this model is 
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To make our approach systematic we will modify this 
model introducing n— index and considering the large— iV 
The usual way to introduce n— index to \4>\ 4 term 



limit, 
is 



(3) 



But it happens that the spectrum of fluctuations around 
the Abrikosov lattice is not positive definite in this case, 
i.e. the Abrikosov lattice is unstabldj. (We will not 
consider the possibilitv-.of condensation into different n- 
components as in RefJj.) To understand the reason of 
this instability let us substitute <j> = <fro + 4>i, where 4>o 
corresponds to the Abrikosov lattice and cf>\ represents 
fluctuations around it, into \<f>\ term. Considering the 
simplest case when fluctuations are small, for term 
we have 



\4>\ = m^l^i + 4>l4>l4>i<t>i + <t>o4>o4>i4>i + 



(4) 



where we wrote only the most interesting, quadratic in <b\ 
part. The spectrum of fluctuations corresponding to (|J) 
was considered by Eilenbergem, and it was found that it 
is positive definite, i.e. the Abrikosov lattice is stable. It 
is important that following to the Eilenberger approach 
one can see that these are the off-diagonal terms in (^) 
which make the lattice stable. But for |0| 4 term defined 
by (|]) in the large— N limit we have 
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so in this case we effectively drop out the off-diagonal 
terms and it leads to the unstable spectrum of fluctua- 
tions. Therefore we suggest the following modification of 
the model: 
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In the large— N limit for \<f>\ term now we have 



where — V M — — acts only in x — y plane with 
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the Abrikosov lattice, and e± = tt^- is the kinetic energy 
corresponding to the motion in the directions perpendic- 
ular to x — y plane. In this Lagrangian cj> is a function 
of coordinates and imaginary time because we want to 
consider the low-temperature quantum regime also. The 
operator \d T \ means that in the Matsubara space this op- 
erator becomes \ui\. For the BCS model the coefficients 
in the Lagrangian (m 
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where eo is the energy of the lowest Landau level, i.e. the 
lowest eigenvalue of m e is the electron mass, H]w 
is the critical field in the mean field approximation, and 
k a = tii /tj_ is the anisotropy coefficient, which is the ratio 
of the in-plane electron hopping and the hopping in the 
direction perpendicular to the x — y planes. 



m<j)04>* n <t>n - 4>l<l>* <i>m<t>m ~ <t>0 00 4>*m 4>*m ) 

which matches the form (^) after the redefinition <p m — » 
i&mi 4>m ~ * — *0 m - It ^ s important that the second term in 
(||) should have the negative sign because otherwise the 
fluctuations around the Abrikosov lattice are unstable 
(from the above simple argument one cannot see that it 
should be negative). So the Lagrangian density of the 
model which we will consider is: 



— + e(p ± ) + \8 T \+ a U n + — - — 
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and the action and the partition function are 
S = [ drd d x£, 
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For simplicity, in the following we will neglect the fluc- 
tuations of the magnetic field, because in the lowest 
Landau level approximation it gives just a renormaliza- 
tion of u termJa'B moreover the high-T c superconductors 
are extremely type-II superconductors therefore even this 
renormalization is not essential. 



III. LARGE TV EQUATIONS. 

The interaction term in the Lagrangian (|^) can be de- 
coupled with the help of a real field p and a complex one 
A 

L = <P* n (£ + |cV| + a + 2pi) cf> n + A>„0„ + AO,* 

B 2 BH N 2 N ... 
+— - — + — p 2 + -A*A, (9) 

07T 47T ZU U 

where £ = — ^— + e(pj_) is the kinetic energy operator. 
From the form of this Lagrangian it is evident that the 
mean field value of p should be imaginary, which is not 
a problem because one can always shift the contour of 
integration over p in the complex plane. Therefore we 
simply redefine p — > —ip, because we will not consider 
the fluctuations of A and p fields. Integrating over TV — 1 
<t> fields we get an effective Lagrangian 

£ = <%(£ + \d r \ +a + 2p)(f> Q + A*0 o ^o + &<P* <P*o 

N n 2 I N A*A I TV In r- 1 I B 2 BH 

where the Green function G satisfies the equation 



£ + M+ a + 2p 2A 

2A* £* + M + a + 2p 



= 6(n-r 2 ) (10) 

Taking the variation with respect to A, p, <j>$ and ne- 
glecting fluctuations of the magnetic field we get the 
large— N equations 
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+ a + 2p(r) Mr) + 2A(r)0*(r) = 0, (12) 



where V± denotes the volume perpendicular to the planes 
with the Abrikosov lattice. 



IV. NORMAL STATE. 

At temperatures higher than the temperature of the 
Abrikosov transition we have cj>o = 0, A = 0, and the 
equations ( p]JT^ ) are reduced to 



1 P= S~ E G(0,0,w,p±), 
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where the limit N — ► oo was taken. The Green function 
in this case is determined by 



(£ + a + 2p + \u\) G(r, 0, u,p±) = 5 2 {r). 



(14) 



Under the approximation when <p field belongs to the low- 
est Landau level we can substitute — tt—D 2 = en, where 
en is the energy of the lowest Landau level and get 
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where e$ + 2p + a = 5 and $n = ■ Also we introduced 
the ultra-violate cutoff e. The advantage of this cutoff 
procedure is that the Green function with the cutoff can 
be written as 
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The phase transition happens when the correlation 
length diverges, i.e. when S — > 0, therefore to find the 
phase transition line we should calculate the integrals in 
( |i~5| ) for the case T 8 in which we have 

4uH /m± \ ^2* 
fi — a* — en = I 1 
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t r(l -d_ L /2) < 5^/ 2 - 1 T, d± < 2, 
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where a* is the "mass" term a renormalized by the quan- 
tum fluctuations 
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£ is a constant 



2 (cothA- ^), 



and we introduced a notation for the dimensions per- 
pendicular to x — y planes d± — d — 2. One can see 
that in high dimensions d± > 4 the correction to the 
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phase transition line is analytical in T, therefore it does 
not change the mean field result essentially. When 
2 < d± < 4 the correction is non-analytical. Substi- 



€q(H — H^*) / H [ c2 ' , for the last case we 
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tuting a* + eo 
have 

H C 2 

where H^* is the mean field upper critical field renor- 
malized by the quantum fluctuations (see (jig)). In the 
physical case d± — 1 one can see that the r.h.s of ( |l7| ) di- 
verges as ^= when 5 
in the form 



0. In this case ( p~7| ) can be written 
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where kq is the Ginzburg number 
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Using (|^) we can estimate kq ~ , where 5 1 is the area 

corresponding to the unit flux $o- By the order of mag- 
nitude S ~ £ 2 , where £ is the coherence length. In case 
when kg is small one can find the qualitative crossover 
line between the Gaussian region (where the u-term is 
not important) and the non-Gaussian region (when it 
becomes important). In the Gaussian regime one can ne- 
glect the r.h.s. in (n9[) getting S « a* + eo- The crossover 
line corresponds to the situation when the r.h.s of jl9| ) 
becomes of the same order with the l.h.s., so that we have 
S - (Vi^T) 2 / 3 or 



H rr — H, 



(0), 
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where H cr is the magnetic field corresponding to the 
crossover from the Gaussian to the non-Gaussian regime. 
The fact that one cannot reach the ordered state going 
from the normal state gives us a hint that the transition 
to the ordered state is not continues. 



V. ORDERED STATE. 



We will solve Eqs.(|10|]l 

when the mean field 



12j) under the approximation 



'o and the fluctuations around it be- 
long to the lowest Landau level. Also, as usually, we will 
consider a lattice with the triangular symmetry. There- 
fore, following to the Eilenberger notations^, we take </>o 
to be proportional to 



(r|0) = (2 ?7 )^, 



(22) 



where 77, £ are components of the vectors which determine 
the unit cell 



r 7 = (1,0), r// = (£,»,) = (i (23) 

In the above formulas we measure all distances in the 
distances between the vertices which is related to the 
magnetic field by 



(24) 



Applying the magnetic translation operator to (|22|) we 
can get a complete basis of functions belonging to the 
lowest Landau level 

0(r|r o ) =e 2 "V^(r + r |0), (25) 

which obey the orthogonality relation 

j d 2 r0*(r|n)0(r|r 2 ) = ^(n ~ r 2 ). (26) 

Restricting our solution to the lowest Landau level 
we artificially decrease the space of solutions of 
Eqs.( [io| , fn]]T2} ), therefore to have a solution we should 
decrease the space of equations projecting them on the 
lowest Landau level. The way to do that becomes evi- 
dent if we consider the Schroedinger equation which cor- 
responds to the Green function equation (nfj) 
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Indeed, to project this equation on the lowest Landau 
level one should multiply it by 0*(r|r o ) and integrate 
over r. Introducing the following notations for the matrix 
elements of p and A 



(28) 
(29) 

(30) 



j d 2 r0*(r|ri)p(r)0(r|r 2 ) = V 2 Pn,r2, 
J d 2 rr(r|rOA(r)r(r|r 2 ) = V 2 A rl . r2 , 
and presenting <fi as 

<Xr) = J d 2 r a{r )cj)(r\ro), 
we get 

(e + e(px) + M + a) a(ri) 

+2 J d 2 r 2 a(r 2 ){pri,v2 + A rl , r2 ) = #a(r x ). 

From the symmetry of this problem we expect that p is 
a periodic function, and A is a quasiperiodic one (that 
means periodic up to the phase), therefore the only ma- 
trix elements which are allowed by the translational sym- 
metry are 



Prl,r2 = Prl^^l - r 2 ), 

A r i, r2 = A rl 5 2 (r! +r 2 ). 



(31) 
(32) 
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Now we can find the eigenfunctions and the energy 
levels of ©: 

«±(r) - -j=^ (e lfcz 0(r|r o ) ± e~ ik ^{v\ - r )) , (33) 
v ± (r) = -±=e^ (ze^0(r|r o ) T ie~ ikz ^{v\ - r„)) , (34) 



E±(r, P± , lo) = e + e±(p±) + a + 2(p r ± |A r |), (35) 
where 



|A r 



One can see that there are two energy branches which 
are denoted by the subscript ±, and there are two eigen- 
functions u and v for each energy level. 

The green function ( |i"o| ) can be expressed through the 
eigenfunctions (^,H) 



G(ri,r 2 ,RL,cj) 
- _J_ v f ■ fl 



^■,,(ri|r )^ s (r 2 |r ) 



(36) 



where the index j = +, — denotes the spectrum branch 
and s — u,v denotes the type of function: 



»*±,u = ( u *± u ±) > 



(37) 
(38) 



The prime under the integral in ( pq ) denotes the integra- 
tion over the half of the unit cell. From the Goldstone 
theorem we expect a singularity in (|3^) at the minimal 
energy, i.e. E- (0, (LO) = 0. This condition is consistent 
with the equation (|lj) if we take 

O = n/WV>(r|0), 

where b is a real positive number. Indeed, projecting jl^ ) 
on the lowest Landau level we get 



e + a + 2( Po - A ) = 0, 



(39) 



which is the same with E_ (0, 0, 0) = 0. 

Now, to get a closed system of equations, we should 
take the matrix elements of the equation ([fl]) .-Using the 
addition theorem of the Eilenberger functionsa. 
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where 



sgnj 



1 for (+) branch 
— 1 for (—) branch 



and we introduced the following functions 



K(r) = 0(r/2|O) 



(r/2|rj/2) 



Q(r) = 1 (0*(r|O)^(O|O) + 0*(r|rj/2)0(O|rj/2) 



(41) 
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From Eq. (|40|) one can see that A r is a linear combina- 
tion of the functions </>*(r|0) and 0*(r|r//2). Therefore, 
presenting A r as 

A r = i (A!r(r|0)^(0|0) + A 2 ^(r|r 7 /2)^(0|r J /2)) , 

(45) 



from (Un) we get two equations on Ai, A 2 
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From the symmetry of the problem we expect that A r 
is not only quasiperiodic (which is evident from ((45|)), 
but it is also quasisymmetric under rotations on tt/3. 
The later property is satisfied only when Ai = A 2 . But 
taking Ai = A 2 we should prove that Eqs.( p6| , p7| ) become 
linear dependent. One can check that it happens indeed, 
multiplying @ by </>(0| ir/) and @ by -<£(0|0), adding 
them, and using that the function 
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D = (V(r|O)0(O|O) + 0*(r|±r,)0(O|±r,) 

x L{v\Q)~m\*i)-kA\*i)mU)) (48) 
transforms under rotation on 7r/3 as 

D -> L> e _i ^ when a; + — ► (x + iy)e i7:/3 . (49) 
Finally we get the following system of equations: 
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(51) 



£;±(r,p_L, lo) = e + a + e(p ± ) + \u\ + 2 (p r ± A|Q(r)|) 



£_ (0,0,0) = 0. 



(52) 



(53) 



VI. SOLUTION OF THE LARGE -N EQUATIONS 
IN HIGH DIMENSIONS. 



In high enough dimensions we can expand Eqs.(j50 51) 
in p, A. This expansion is possible when either d± > 4, 
T 7^ or d± > 2, T = 0. In these cases one can write 
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From Eqs.@||) we get 

- 4/3A J d 2 r Q*(r)Q(r) - (bu - A)Q(0), 
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(55) 



2 / d 2 r {a - 0(e o + a + 2p(r ))) K(r Q - r) 

= p r -buK{r). (56) 
The first equation ( |55"| ) can be solved immediately giving 
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The second equation (|56j) can be solved by the Fourier 
transformation which wc define as 



K(v) 



k 



K(k), 



(58) 



where due to the periodicity of K(r) the vector k should 
have the following discrete values 



k = 

,k 2 



kini + k 2 n 2 , 
= 2.(1,-^), 



(59) 
(60) 



where ni, n 2 are integers. The Fourier transformation of 
K(r) can be done analytically 



(61) 



Defining the Fourier transformation for p T by the same 
rule, from Eq.(Bq) we get 



P(k) 



4(q - /3(e + a))r) 5 K0 + buK(k) 
1 + i/3r]K(k) 
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Substitution of ( |57| , |62|) to ([33]) gives the equation for the 
condensate density b 
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Using this equation we can present the spectrum as 
v K(k) e lkr |Q(r)| % 



E±(t) = bu 
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(r = Q). 



(64) 



One can check that this spectrum is positive definite 
which means that our solution is stable. The transition 
point can be found from (|6^) taking b = 



eo + a c = -877a, 



(65) 



where a c means the critical value for a. The same result 
can be obtained from (JT^) . We expect to have a nonzero 
b in the ordered state when a < a c . But Eq.(|63[) has posi- 
tive solutions for b (and b was chosen to be positive) only 
when the expression in the parenthesis on the l.h.s. of 
d63| ) is negative. Numerical calculation of this expression 
gives that it is negative only when 



/3</3 c 



(66) 



where f3 c — 0.112. When the condition (|66| ) is satisfied 
the condensate density b can have any small values as a 
function of a. Therefore in that case the phase transi- 
tion is of the second order. Note that the expansion (|54|) 
is valid only when A, p, b <C eo- When the phase transi- 
tion is of the second order this condition can be always 
satisfied if we close enough to the phase transition line. 
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Therefore the absence of the solution when (3 > (3 C in fact 
means that the transition is not continues. 

Let us now summarize the results of this section where 
we were interested in the cases d > 4, T ^ and 
d > 2, T = : When the interaction constant is not 
too strong (condition (fi^)) the phase transition is of the 
second order. If the condition (|6^) is not satisfied then 
the transition between the normal to the ordered states 
cannot be continues. 



VII. PHYSICAL CASE 

Let us consider the physical case d = 3 (d± = 1). Tak- 
ing the integral over p± and summing over to in ( j50|]5l| ) 

we get 



t / dfr 



/T(/(e + (r)/T)-/(e_(r)/T) 

VMO)|lQ(r)l 



(bu - A)Q(O) (67) 



t / d 2 r Q 



>/T(/(e + (ro)/T)-/(e_(ro)/T) 



— (VMro) - yje-(r ) 



K{r - r) 



where 



, . 4V3l 
p r — buKyr) 



e±(r) = e + a + 2(p r ±A|Q(r)|), 



-Xx 



coth 7rA 



and 



2t] 2 1 2 V2tt 



(68) 



(69) 



(70) 



(71) 



It is convenient to introduce the following dimensionless 
variables 



T = e K 2 G T', A = eoK%A', 



(72) 



fc / d 2 r 



<T 
2 



7 (/( e V(r)/T')-/( £ '_(r)/T')) 



|Q(r)| 



(b' - A')Q(O), (75) 



jfe / <i 2 r 



'(/( e ' + (ro)/T')-/(e'_(r )/T')) 
K(r - r) 



ire 



(76) 



where fc = l/(2v37r) is a numerical constant. Note that 
the third term on the r.h.s. of Eq.(|76]) can be absorbed 
into renormalization of p' and a' 



7re 



-t 



a — > a 



8V3e 
— — < 

7re 



This corresponds to the renormalization of the external 
magnetic field by the quantum fluctuations as in (111). 
Eq.([76j) is an integral equation with the kernel 



K(r) = e -^|(™?+™i™2+™2) e -27ri(kini+k2n2)r 



where It! , k 2 were defined in (|60|) . To simplify the prob- 
lem we will take an approximate expression for this kernel 
leaving only zeroth and first harmonics: 



K(t) -2 + 7 



n 2y „ / y \ 

cos 2-n—= + cos Z7T I x H = 

V3 V VV 

■ cos 2tt [ x -= ) 

V3/ 



(78) 



where 7 = 4e~ 27r / v/ ^ w 0.1063. It seems to be a good 
approximation because the contribution from higher har- 
monics to (77) decreases exponentially with the harmonic 
order. For example the contribution from the second har- 
monic is less then 0.1% of the first one. Presenting p' r in 
the form 



p = e K 2 G p' , bu = e K G b', 



2 u 



(73) 



(74) 



Pr = P0 + 7 Pi 



2y ( y \ 
cos 2n—= + cos 2n \ x H = 

V3 V VsJ 



■ cos 2n ( x -= ) 



(79) 



where kg is the Ginzburg number (|20|). In these new 
variables the equations do not contain kq explicitly: 



we get the following system of equations 
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2k / d 2 r 



\f{e' + {v)/T') + f{e'_{v)/T') 
= Po- 26' 



( x /eV(r) + ,/e'_(r) 



fc / d 2 r 



IT 



'(/(e' + (r)/T')+/(e'_(r)/T')) 



'(/( e ' + (r)/T')-/(e'_(r)/T') 



(80) 



cos2tt-^= = pi - b', (81) 



|Q(r)| = (&' - A')Q(O), 

(82) 




: (r) = 2(p l7 afoy) + A'(Q(0) ± |Q(r)|)) (83) 



p = A'Q(O) - 3 Pl7 - o'/2 (84) 



where 



s(x, y) — cos 2n—^= + cos 2tt ( x + -K=] 

+ cos2tt (x - -^=J - 3. (85) 

Eq.(||) follows from Eq.(|5^). Note that the spec- 
trum contains only A' and pi, therefore combining 
Eqs.( |31| , |B2| ) to get rid of b' on the right hand sides we get 
a closed equation on A' and p\ . Therefore taking a given 
A' we can solve this equation for p\. Then, knowing A' 
and p\ we can find all the other parameters a', po- In 
Figs. 1-3 we present the graphs of pi, 6' ,a' as functions 
of A' for different temperatures. Note that for large A' 
the fluctuation contribution is small and one should have 
the mean field results: 



b' 



Pi 



A', 



-2Q(0)6', A'-voo. 



(86) 



One can see that the graphs on Figs. 1-3 begin to ap- 
proach this asymptotics. 

Let us start the analysis of these graphs from the low- 
temperature case (see Fig. 1): One can see that a' first 
increases as A' decreases (as one should expect from the 
mean field theory), but then it decreases. The extremum 
point corresponds to the first order phase transition. 

At high temperatures the behavior is different (see Fig. 
3): Going from large A' we see that in this case first the 
condensate density b' becomes zero, and then there is 
an extremum in a'. In our approach it was chosen that 
b' > 0, therefore the point b' = corresponds to the first 
order phase transition. We call this transition also first 



FIG. 1. Solution of Eqs.(g0[j84|) at temperature T' = 0.1. 
The extremum in a 1 corresponds to the first order phase tran- 
sition. Note that the condensate density b' is finite at this 
point. 



order one because it is still discontinues in p and A. So, 
the transition is of the first order for all temperatures, but 
at low temperatures it corresponds to the extremum of 
a', while at high temperatures it corresponds to depleting 
of the condensate density to zero. This change in the 
kind of the first order phase transition takes place at 
T = T* w 2.6. 

The phase transition curve which was obtained from 
the above criteria is 



a> c (T') = h(T>), 



(87) 



where h is plotted on Figs. 4,5. The asymptotic behavior 
of the function h at large arguments is 

h{x) w -2.92 x 2/3 + 2.31 x 1/2 , x > 1, (88) 

with an accuracy of leading and sub-leading terms. At 
low temperatures the function h(T') is analytical in T'. 
In spite of the fact that there is a change in the kind of 
the transition at T* , this curve does not have a break 
at this point. The phase transition line in the original 
notations is 



T 



(89) 



where a* is a renormalized by quantum fluctuations (|1^ 
One can see that this curve scales with the Ginzburg 
parameter kg- 
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A" 



A" 



FIG. 2. Solution of Eqs.(|8Cf|84|) at temperature T" = 2.0 
The condensate density b' is close to zero at the point of 
extremum of a' . 



FIG. 3. Solution of Eqs.(gfjj84|) at temperature T" = 20. 
The first order phase transition corresponds to the point 
where b' = 0. 



According to the form of ( |89| ) one can define too 
regimes: the low-temperature, quantum regime 



T < € Kq, 

and high-temperature, classical one 
T > £ 4- 



(90) 



(91) 



Note that the quantum region extends as kg increases. 
Asymptotically in the classical regime from rt§9) one gets 



a* +e = -2.92 e 



which agrees with the estimation (pl|). 



2/3 



(92) 



VIII. LOW-ENERGY SPECTRUM. 

In this section we will show that the low-energy spec- 
trum of fluctuations in our model is different from one 
which follows from the Eilenberger theory. In our model 
the Eilenberger result can be obtained if one neglects the 
terms on the l.h.s. of Eqs. (|5C|j5ll ). These terms are small 
if we are not too close to the phase transition line, never- 
theless they are always not equal to zero, and as we will 
show they change the low-energy spectrum considerably. 
If we neglect the mentioned terms in Eqs. (|5^^l|) then we 
have 



A = bu, 
p r = buK(r), 



(93) 
(94) 



and for the lower energy branch we get 

E_(r) = 2bu(K(r) ~ \Q(r)\) - 2bu(K(0) - |Q(0)|), (95) 

where we used that -EL(O) = 0. It happens that if we 
expand this expression in powers of r, then the terms 
quadratic in r cancel each other and the expansion starts 
from r 4 term 



E_(r) 



r 4 . 



(96) 



This result leads to divergences, for example the fluctu- 
ation contribution to the density in case d = 3 ,T =/= is 
logarithmically divergent 



d 2 r dp 
p 2 + E- (r) 



d 2 r 



(97) 



Also the fluctuation ccpiection to the conductivity has a 
logarithmic singularity^. The fact that there are infra- 
red divergences means that one needs a more careful anal- 
ysis of the infra-red behavior of the model. The similar 
situation happens in the two-dimensional Bose gas at non 
zero temperatures. The Bogolubov approximation leads 
to the divergent fluctuation contribution to the density, 
but the careful analysis of the infrared behavior based 
on the effective low-energy functional leads to the theory 
without divergences. To the best of our knowledge the 
asymptotic behavior of the model under consideration 
was not found yet. 

In our large- model this problem does not arise be- 
cause due to the fluctuation contribution the fine tuning 







FIG. 4. First order phase transition line at low tempera- 
tures. 



FIG. 5. First order phase transition line. 



between the particle-hole and particle-particle parts of 
the spectrum (p r and A|Q(r|) does not happen, and one 
expects that the leading term in the low-energy spectrum 
is r z 



E-(r) 



(98) 



In case of high dimensions one can see this explicitly 
from Eq.(64). And in general case, one can see that the 
particle-particle and particle-hole parts of the spectrum 
are affected by the fluctuation terms in a different way 
and therefore we do not expect any fine tuning between 
these terms. 

We think that the result E_ (r) ~ r 2 is specific for the 
large— A model, and the situation in the real model is 
much more complicated. Nevertheless, we think that the 
large— A model is a reasonable model for the description 
of the phase transition, because the infra-red properties 
seem to be irrelevant for the phase transition. Indeed, 
usually the infra-red divergences are absent in the per- 
turbation expansion for the physical quantities like free 
energy, density, etc. And it is enough to know the free 
energy to determine the kind of the phase transition and 
to find the phase transition line. 



IX. DISCUSSION AND CONCLUSIONS. 

We considered the effect of order parameter fluctua- 
tions on the transition between normal and mixed super- 
conducting states in pure superconductors. Our start- 
ing point was an effective functional of GL type. We 



showed that the coefficients in this functional are finite 
(i.e. this functional exists) in the quasi-two-dimensional 
situation, when the applied magnetic field is parallel to 
the low-conducting direction. This case is interesting 
from the point of view of high— T c superconductors, be- 
cause they have a quasi-two-dimensional band structure. 
We considered this functional in the large— A" limit. One 
should be careful introducing the n-index into the La- 
grangian, because the symmetry between the particle- 
hole and particle-particle channels is very important for 
this problem. Indeed, introducing the n-index in the 
usual way 



(99) 



in the large— A limit one effectively drops out the 
particle-particle channel, and it leads to the model with 
an unstable spectrum of fluctuations. Therefore we in- 
troduce the n-index in the following way: 



4>* "(jjcj)* "(j) ->■ 2(j)* n (j) n (j)* m (j) m - 4>* n cj)* n (j) m (j) m . 



(100) 



The coefficients in the above formula can be found from 
the following consideration: The effect of fluctuations can 
be formally suppressed reducing the Ginzburg number. 
And in this limiting case the Eilenberger theory becomes 
exact. We want our large— A model to be as close to the 
real model as possible, and therefore in the limiting case 
kq — > we should have the Eilenberger answer for the 
spectrum. Th is requirement uniquely defines the coeffi- 
cients in (100). 

To simplify the large— A equations we used the lowest 
Landau level approximation which is valid when the order 



10 



parameter is much smaller than eo, i.e. when we are not 
too far from the phase transition line. These large— N 
equations can be easily solved in case of high dimensions: 
either d± > 4, T ^ or d± > 2, T = 0. In these cases 
the transition was found to be of the second order if the 
interaction constant is not too large. It is interesting 
to draw a parallel between our solution of the large— N 
model and the renormalization group, approach to the 
quantum critical phenomena problemsliireL3. In our case 
the dynamical exponent z = 2. Note that the magnetic 
field "eats" two dimensiopSp-fcherefore the straightforward 
application of the resultstirll3 gives that the upper critical 
dimension at zero temperature is dj_ = 4 — z = 2, which 
agrees with our approach. The phase transition line in 
the large— N limit was found to be 



(())* 



H C 2 - H c2 

r(0) 



-T 



d±/2 



2<d±<4, 



(101) 



where H c2 is the mean field upper critical field renor- 
malized by the quantum fluctuations. In fact this result 
is more general than that of the large— N limit. Indeed, 
the difference between our model and the standard GL 
model (which was considered in Rcniil) ariseSrOnly when 
one considers the renormalization of u termd. But in 
the case under consideration the u term is irrelevant and 
therefor e one should get the usual answer. Therefore the 
answer (101) should hold in case N = 1 too. 

In case of physical dimensionality (d± = 1) the model 
gives the first order phase transition. The fact that at 
finite temperatures the phase transition is of the first 
order looks natural because the fluctuation contribution 
diverges as one reaches the phase transition line from the 
normal state. (For example the first order correction to 
the "mass" term diverges as 1/VS, see (fl7|).) This situ- 
ation is similar to one which happens in the model stud- 
ied by Brazovskiy in RefJj, where the fluctuations drive 
the phase transition to the first order one. Therefore we 
think that in the real model at finite temperatures the 
transition is of the first order too. At zero temperature 
the large— N model also gives the first order phase tran- 
sition, but it is not clear whether in the real model the 
transition should be necessarily of the first order at zero 
temperature. 

The phase transition line which follows from our model 

is 



a* + e = € Kq h 



T 



e Q K 2 G 



where h is plotted on Figs. 4,5, and 



co 



(0)* 



H c2 - H c2 



H. 



(0) 



(102) 



(103) 



Note that we considered only the low-temperature part 
of the phase diagram T <C eo, so that the result (103) 
may be applied only in this case. According to the form 
of (103) one can define to regimes: T -C eo^g corre- 
sponding to the quantum fluctuations, and T 3> £qKq 



corresponding to the classical ones. Note that increase 
of the Ginzburg number makes the problem more quan- 
tum. For example if kg ~ 1, then one cannot reach the 
classical region because our theory works when T -C eo- 
Asymptotically in the classical region the phase transi- 
tion line is 



a* +e Q = -2.92 e 
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(104) 



Qualitatively the phase transition line looks similar to 
the experimental datatj on overdoped high— T c materials: 
The upper critical field significantly increases as temper- 
ature decreases showing a nonanalytical dependence. In 
our theory, the curvature of the phase transition line is 
negative in the classical regime, but at low temperatures 
it becomes positive (see Fig. 4). Note that the Ginzburg 
number in this problem is proportional to the anisotropy 
K G ~ "Tci that enhances the fluctuation contribution. 

Indeed the resistive phase transition is broad in this ma- 
terials, that supports that the fluctuation contribution is 
large. 

Finally, to avoid confusion, we note that in the high 
temperature superconductors, a line in the H-T plane re- 
ferred to as the irreversibility line, is usually interpreted 
in terms of the melting of the vortex lattice. To ad- 
dress these experiments our work should be generalized 
to incorporate the effects of disorder which are beyond 
the scope of our paper. However, low enough disorder 
should not affect the melting line. Therefore in that case 
the transition line which was found in the paper can be 
considered as the irreversibility line. 
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